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The multiconfiguration Dirac-Hartree-Fock theory (MCDHF) has been employed to calculate the 
electric dipole moment of the 7s6d 3 D2 state of radium induced by the nuclear Schiff moment. The 
results are dominated by valence and core-valence electron correlation effects. We show that the 
correlation effects can be evaluated in a converged series of multiconfiguration expansions. 
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I. INTRODUCTION 



A non-zero permanent electric dipole moment (EDM) 
of an atom, molecule, or any other composite or ele- 
mentary particle is one of the possible manifestations of 
parity (P) and time reversal (T) symmetry violations. 
In the absence of any external electromagnetic field an 
atom can have a permanent EDM either due to an in- 
trinsic EDM of one of its constituent particles or due to 
P- and T-violating (P-odd and T-odd) interactions be- 
tween these particles 0, . When compared with the in- 
trinsic EDM of the constituent particles, the net induced 
EDM of an atom or molecule is often expected to be 
larger by several orders of magnitude due to various nu- 
clear and atomic enhancement mechanisms. Therefore, 
atoms and molecules arc considered to be very attrac- 
tive for carrying out EDM experiments and in the search 
for 'new physics' beyond the standard model, since, in 
the latter case, the induced EDM is greatly suppressed, 
when compared to the anticipated values from the 'new' 
theories [3|. In atomic physics, in particular, the experi- 
mental search for a permanent EDM is gaining momen- 
tum due to recent advancements in trapping free neutral 
atoms [1, H, Hj] , including various radioactive species 



During the last decade several atoms were considered 
as candidates for such experiments [1, Q. These in- 
volved (i) diamagnetic atoms (i.e., total angular mo- 
mentum J — 0) in their respective ground states, (ii) 
the alkalis, and (iii) atoms with a single p electron 
outside closed shells, which were investigated in laser 
traps for their prospectives to perform EDM experi- 
ments. Presently, radium appears to be the most promis- 
ing candidate, and experiments on this element are under 
way at the Argonne National Laboratory [7, 1C . [Til as 
well as the Kernfysisch Versnellcr Instituut 12L 13l. Il4|. 
The main advantages of radium lay in (i) large nuclear 
charge Z , (ii) simple electronic structure and closed-shell 
[Kr]4d 10 4/ 14 5s 2 5p 6 5d 10 6s 2 6p 6 7s 2 % ground state, (iii) 
octupole deformations of the radium nuclei for several 
isotopes [H, IH| > as well as in (iv) coincidental proxim- 
ity of two atomic levels of opposite parity, 7s7p 3 Pi and 
7s6d 3 D2, which are separated by a very small energy 
interval 5.41 cm . In particular, the latter two advan- 
tages give rise to a relatively large enhancement factor, 
which is one of the largest among the atoms considered 
so far P3, [3. 

However, extraction of fundamental P- and T-violating 
parameters or coupling constants from experimentally 
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measured atomic EDM requires atomic form factors 
which can be provided only by an ah initio atomic the- 
ory. Several of these form factors have been previously 
calculated by the group of Flambaum 0, QjJ ■ In practice, 
there are essentially four different form factors, related to 
four mechanisms which can induce an atomic EDM. An 
atom can acquire a permanent EDM due to P- and T- 
odd clectron-nucleon interactions, or due to the electro- 
magnetic interaction of atomic electrons with nuclear P- 
and T-odd moments, of which the leading ones are: the 
Schiff moment, the magnetic quadrupole moment, and 
the electric octupole moment, respectively. The latter 
two moments may exist only in nuclei with spins larger 
than / = 1/2, while the Schiff moment may exist also 
in isotopes with nuclear spin / = 1/2. Nuclei with spins 
larger than 1=1/2 produce electric quadrupole shifts 
which arc difficult to account for in an EDM measure- 
ment. Therefore an EDM induced by the Schiff moment 
in 1=1/2 isotopes seems to be the property of choice 
among (most of) the experimenters [13l [l9l . [20j . The 
nuclear Schiff moment is a P-odd and T-odd (elcctric- 
dipolc) moment that occurs due to P- and T-violating 
interactions at the nuclear scale. The Schiff moment 
mixes atomic states of opposite parity and may induce 
static EDM in atoms if magnetic and finite-size effects are 
taken into account in the electron-nucleus interaction Q . 

In the present paper, we present calculations for the 
atomic EDM in radium, as induced by the Schiff moment. 
The atomic wave functions were obtained within the 
framework of the multiconfiguration Dirac-Hartree-Fock 
(MCDHF) theory. The wave functions were separately 
optimized for the 7s6d 3 i?2 and 7s7p 3 Pi states (simi- 
lar calculations were carried out recently for the scalar- 
pseudoscalar contribution to the EDM in cesium [2~0]). 
The main purpose of this paper is to provide a systematic 
evaluation of the effects of electron correlation on the cal- 
culated EDM of radium. We demonstrate the saturation 
of the core- valence correlations, the dominant electron 
correlation effect beyond the Dirac-Fock approximation. 



II. THEORY 

The Hartree-Fock and Dirac-Hartree-Fock theories, 
based either on the finite-grid, basis-set, or some other 
numerical methods provide a natural point of depar- 
ture in describing the electronic structure of atoms and 
molecules. For medium and heavy elements, these 
methods are often combined with Breit-Pauli or Dirac- 
Coulomb-Brcit Hamiltonians in order to account for rel- 
ativistic and retardation effects on the wave functions 
and the level structure of complex atoms. However, the 
main hindrance in applying modern computational tech- 
niques arises from the electron-electron correlation, i.e., 
the residual interaction among the electrons beyond the 
atomic mean field, and this is especially true for sys- 
tems with many electrons. In neutral or nearly-neutral 
systems, missing electron-correlation effects are indeed 



often the main reason for the discrepancies between the 
observed and calculated properties of atoms. 

Today, many-body perturbation theory (MBPT) [22], 
HDjIULHII and various variational methods, often referred 
to as the multiconfiguration Hartree-Fock (MCHF) the- 
ory (or its relativistic counterpart — the Dirac-Hartree- 
Fock theory [26|, [2?], Hfl ) are the two dominant pillars in 
performing atomic structure calculations. These meth- 
ods are designed to evaluate the electron correlation ef- 
fects in a systematic manner. For MBPT and most re- 
lated methods, a nearby closed-shell configuration of the 
atom or ion is typically a convenient starting point. Dur- 
ing the past decades, therefore, MBPT techniques were 
mainly applied to systems with either no or just a single 
electron outside closed shells. Systems with several elec- 
trons in open shells arc, in contrast, much more difficult 
to deal with if benefit is to be taken from the theory of 
angular momentum and spherical tensor operators, i.e., 
by using a restricted representation of the one-electron 
orbitals. Difficulties occur then due to the large depar- 
ture from a closed-shell V N potential and the rapid in- 
crease in the complexity of all perturbation expansions 
with any additional electron outside closed shells. Al- 
though in the variational multiconfiguration methods the 
algebraic complexity also depends on the shell structure 
of the atom, these methods can be applied more easily 
to systems with an arbitrary number of electrons outside 
closed shells. Apart from the number of open shells, the 
accuracy of multiconfiguration Dirac-Hartree-Fock calcu- 
lations depends crucially also on the occupation of the 
valence shell(s) but this occurs rather indirectly, through 
the limiting number of configuration state functions that 
can be included in a particular wave function expan- 
sion. In this sense, the limitations are less conceptual 
but arise from the available computer resources. These 
limitations are typically related to the structure of the 
valence shell(s), i.e., to the angular properties of the va- 
lence electrons and their couplings. 

A. MCDHF theory 

We used a slightly modified version of the General Rel- 
ativistic Atomic Structure Package (GRASP) [28| to gen- 
erate the electronic wave functions. In the multiconfigu- 
ration Dirac-Hartree-Fock method, the wave function for 
a particular atomic state ^('jPJMj) is obtained as a lin- 
ear combination of configuration state functions (CSFs) 
which are eigenfunctions of the parity P and the total 
angular momentum operators J 2 and J z , 

NCF 

* {jPJM.j ) = ^] C r $ (j r PJMj ) . ( 1 ) 

r 

In the present computations, the wave functions were 
separately generated for the 7s7p 3 P± and 7s6d 3 L>2 
states of radium. Each wave function was obtained as 
self-consistent solution of the Dirac-Hartree-Fock equa- 
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nuclear radius [fm] 



FIG. 1: Nuclear charge density p [a.u.] as function of the nu- 
clear radius r [fm], using a two-parameter Fermi distribution 
for the 2 g|Ra nucleus. 



tions [26j by using a systematically increased multiconfig- 
uration bases (of size NCF) of symmetry-adapted config- 
uration state functions &(-f r PJA4j). Configuration mix- 
ing coefficients c r were obtained through the diagonali- 
sation of the Dirac-Coulomb Hamiltonian 

N 

H DC = ]T [coy • Vi + (ft - l)c 2 + V(r 3 )] + £ l/r jk 

j=l j>k 

(2) 

where V(r) is the monopole part of the electron-nucleus 
interaction. A more detailed description of the the- 
ory [27], Us} and method of calculation [3(], HH can be 
found elsewhere. 



In the GRASP code [28| all electronic (radial) orbitals 
are represented on a numerical grid which increases expo- 
nentially in order to ensure an accurate representation of 
the atomic wave functions near the nucleus. The grid 
is generated from the formula r.; = r D exp((i — l)/i), 
with r a = 2.0 x 10~ 8 x a , h = 7.0 x 10~ 3 , and 
i = 1, . . . , 4000. With the above parameters, there were 
1255 grid points within the 'half-chargc-density' nuclear 
radius 7' nuc = 6.85 fm, while the numerical represen- 
tation of the full extent of all radial (core and valence) 
electronic orbitals required more than 3000 points. With 
the above choice of the radial grid, all necessary onc- 
and two-particle matrix elements can be calculated with 
a (relative) accuracy of the order of ~ 10 -8 or better. 
Note however that larger uncertainties may arise from 
the radial matrix elements in Eq. (|14[) . due to the ap- 
proximate nature of the nuclear charge distribution in 
Eq. ([3]) and due to deviations from radial symmetry of 
those isotopes, for which nuclear deformations are signif- 
icant (for model dependence see [HI). 



C. Atomic EDM 

Neglecting the contributions from the off-diagonal hy- 
pcrfme interaction the coupled wave function of the to- 
tal system 'electrons + nucleus' is given by the Clebsch- 
Gordan expansion [35| 

* (jvPJIFMf) = 

(JIMjMz\JIFM F ) ^(jPJMj) V(uIMi) (4) 

MjMi 



B. Nuclear shape 

All sclf-consistcnt-ficld calculations were done with the 
nucleus modelled as a spherical ball by applying a two- 
parameter Fermi distribution 

P(r) = — ff- rj-r (3) 

1 + exp[(r — c)/a\ 

in order to approximate the radial dependence of the 
nuclear charge density p(r). The parameter po is de- 
rived from the normalization condition J p(r)d 3 r = Ze. 
Figure [T] shows the nuclear charge density p(r) inside 
the 2 ||Ra nucleus, calculated with the parameters c = 
6.85 fm and a = 0.523 fm (see [H, H3 for details). For 
the other two isotopes the 'half-charge-density' parame- 
ter c was set to c = 6.83 fm ( 2 g|Ra) and c = 6.73 fm 
( 2 ggRa), with the 'nuclear-skin-thickness' a value un- 
changed. These nuclear charge distributions were also 
used in the subsequent calculations of the Schiff moment 
expectation values (Eq. [T4")) . 



where ^ (i/IMj) represents the ground state of the nu- 
cleus, and where the standard notation is used for the 
Clebsch-Gordan coefficients. For high-Z elements with 
closed levels of opposite parity, such as radium in the 
3 Pi and 3 D2 levels, one of the most important parity (P) 
and time reversal symmetry (T) violating interactions is 
caused by the nuclear Schiff moment S, which gives rise 
to the electron-nucleus interaction: 

jv 

Hsu = 4tt £ (S-V,-) p( rj ) . (5) 
3=1 

In this Hamiltonian, p (r) is the normalized to unity nu- 
clear density function from Eq. ((3]), and the Schiff mo- 
ment S is directed along the nuclear spin I. S = SI /I. 
The interaction in Eq. (J5]) mixes states of different parity 
and may also induce a static electric dipolc moment of 
the atom. Since the Schiff moment interaction is quite 
weak, we can express the wave function of the (mixed- 
parity) hyperfine state \F,M F ) of the level J p as pll ]: 
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* ( 1V JIFM F ) = a* (-yvPJIFMp) + ^ h * ( ai v{-P)JJFM F ) (6) 

i=i 

where the coefficient a of the given hyperfine state can be set to 1. The expansion coefficients of the other (hyperfine) 
states of opposite parity can be perturbatively approximated by 

(* (aiv(-P)Ji IFM F ) \H S m\^> (jvPJIFMf)) 

h. = J L (7) 

E(jPJ) - E{oi(-P)Jd ' 

The mixed-parity wave function in Eq. (J6]) for the hyperfine state \F, Mp) of a particular atomic level 2S+1 Lj induces 
a static EDM of an atom: 

m 

D A = (^vJIFM F )\D z \^ (^vJIFMf)) = 2^ bi(^ ^vPJIFM F )\D z \^ (onv{-P)JiIFM F )\ (8) 

where D z denotes the z projection of the electric-dipole moment operator. For this electric-dipole operator, the matrix 
element between (hyperfine) states of different parity can be expressed as: 



* {'yvPJIFMp) \D Z \^ {a t v{-P)JJFM F ) 
^y+J+F+i {2F + 1) V 2JTl {^_ F Mp I M F ){ J FJi {}[*(7^)ll^ 1 ||*(« i (-^)^)] (9) 

while the matrix element of the (scalar) electron-nucleus interaction in Eq. ([5]), induced by the Schiff moment, is 
written as: 

{ 1V PJIFM F ) \H S m\* {aiv(-P)Ji IFM F j) = 

(-1)*+'+™ V2J + T ^Z±E ){ j, ' 1 | 47 r5 [*{ lP J)\\^ p{r)\\*{a t {-P)J t )] (10) 

I 



and is independent of Mp- In the following, we shall 
refer to the last term on the right-hand side of Eq. (flQ|) 
as the reduced matrix element of the Schiff operator for 
the two (fine-structure) levels J and Ji of different parity, 
and shall assume Mp = F, in line with optical pumping 



schemes of hyperfine levels with circularly polarized light. 
For the \F,Mf) hyperfine state of the 3 D 2 level, the 
static EDM in Eq. ©, induced by the P-odd and T-odd 
nuclear Schiff moment, becomes 



D A { 3 D 2 , FM F ) = 



(* { 3 D 2 , FMp) \b z \^ ( 3 P U FM F )) (* ( 3 Pi, FMp) \H SM \* ( 3 D 2: FM F 
2 E( 3 D 2 ) - EjtPO 



(11) 



if the summation over the intermediate states is restricted 
to the nearby 3 P\ level. This assumption is justified by 
the size of the energy denominator which is 500 times 
smaller than the second smallest. Moreover, the electron 
nucleus interaction in equation (O is scalar, therefore 
only one of the hyperfine states may occur. 



When the multiconfiguration expansion from Eq. ((T|) 
is employed for the electronic part of the (total) wave 
functions, the reduced matrix elements of general tensor 
operator can be decomposed into (reduced) matrix 
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elements between configuration state functions 
[*( 7 PJ)||f fe ||*(a(-P)J i ) = 

J2 crc s h { lr PJ) \\T k \\<S> (7 s (-P)Ji)] (12) 



and those, in turn, into a sum of single-particle matrix 
elements 



and 



$(7,PJ)||r fc ||$(%(-^)^) 



J2 d *b( rS ) [ n aKa\\t h 



\n b K b \ 



(13) 



In the latter expansion, the d^ b (rs) are known as 'angular 
coefficients' that arise from using Racah's algebra in the 
decomposition of the many-electron matrix elements (27l . 
l36j . The single-particle reduced matrix elements in the 
expansion (|13j) can be factorized into reduced angular 
matrix elements and radial integrals which, for the Schiff 
moment interaction, read 



n-aKaWV 1 p(r) \\n b K b 



[KaWC^Kb] I (P a P b + Q a Q b ) ^ dr , (14) 



dr 



while for the electric-dipole moment operator (fc = 1) 
it is 



na,Ka\\d 1 \\n b Kb = 



-[/tallC 1 !!^] / (PaPb + QaQb) r dr. (15) 
Jo 

For the calculations of the matrix elements we extended 
the GRASP [28[ and mdfgme [37| relativistic atomic 
structure packages. The extension, presented in this 
work, includes programs for both Schiff moment inter- 
action and electric-dipole moment matrix elements. Ex- 
perimental energy differences were used in the calcula- 
tions of all expectation values. The energy values for the 
two levels of interest (-EV S 6d sd 2 = 13993.97cm -1 and 
^7s7p 3 Fi = 13999. i&cm^ 1 ) were taken from the tables 
of Moore [38[. The nuclear spin and magnetic moment 
data were taken from the tables of Raghavan [3{| . 



D. Handling non-orthogonalities 

The electronic wave functions were optimised sepa- 
rately for the two levels of interest. All expectation values 
were evaluated with the biorthogonal technique devel- 
oped by Mahnquist [40| . For two atomic state functions 



*(7PJ) = ^c r $( 7r P.7) 



*(a(-P) J,) = £><&( 7s (-P)Ji; 



the reduction of a general matrix element 
y{"fPJ)\\f k \\y{a{-P)J t ) 



(17) 



(18) 



(16) 



into a sum of radial integral x angular coefficient terms 
is based on tensor algebra techniques. In the decomposi- 
tion p3p , it is usually assumed that the (many-electron) 
configuration states on both sides of the matrix clement 
are built from a common set of spin-orbitals. This is 
a very severe restriction since a high-quality wave func- 
tion demands orbitals optimized for the specific electronic 
state. Instead of the standard decomposition (fT3"|) based 
on tensor algebra techniques, Malmquist has shown 
that for very general expansions, where the two atomic 
states are described by different orbital sets, it is possi- 
ble to transform the wave function representations of the 
two states in such a way that standard techniques can be 
used for the reduction of the matrix elements in the new 
representation. This procedure has been implemented in 
the modules that compute Schiff moments and it can be 
summarized as follows 

1. Perform MCDHF or CI calculations for the two 
states where the orbital sets of the two wave func- 
tions are not required to be identical. 

2. Change the wave function representations by trans- 
forming the two orbital sets to a biorthogonal basis. 
This is followed by a counter-transformation of the 
expansion coefficients c r and c s so as to leave the 
resultant wave functions invariant. 

3. Calculate the matrix elements with the trans- 
formed wave functions for which now standard 
techniques can be used ■ 

The transformation of wave functions is very fast, since 
it relies only on angular coefficients for a one-electron 
operator of rank zero which appears in the evaluation of 
the kinetic energy term in the MCDHF or CI step. The 
details of the transformations arc discussed in 14111. 



III. METHOD OF CALCULATION 

To generate the atomic states of interest, the method 
described as systematic expansion of configuration set 
[30l . l3lT ] has been employed, in which symmetry-adapted 
CSF of a given parity and total angular momentum 
are generated by substitutions from reference configura- 
tions to an active set of orbitals. The active set should 
hereby comprise (at least) all valence shells, several 'near- 
valence' core shells, as well as certain number of virtual 
shells. The active set and multiconfiguration expansions 
are then systematically increased until the expectation 
value (of interest) is converged. In practice, we divided 
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the computations into three phases to generate (i) the 
spectroscopic orbitals, (ii) the virtual orbitals, and (hi) 
to perform large configuration interaction calculations, 
once the set of orbitals is fixed. 



A. Orbital set 

In the first phase of the computations, all spectro- 
scopic orbitals required to form a reference wave func- 
tion were obtained with a minimal configuration expan- 
sion, with full relaxation. These orbitals were determined 
from a symmetry-adapted Dirac-Fock calculation with 
only those configurations which arise in j-j coupling for 
a particular state of interest. The spectroscopic orbitals 
were kept frozen in all later steps. 

In the second phase the virtual orbitals were gener- 
ated in five consecutive steps. At each step the virtual 
set has been extended by one layer of virtual orbitals. A 
layer is defined as a set of virtual orbitals with different 
angular symmetries. In the present paper five layers of 
virtual orbitals were generated, each layer comprising or- 
bitals with symmetries s,p,d,f,g, and h. At each step the 
configuration expansions were limited to single and dou- 
ble substitutions from valence shells to all new orbitals 
as well as to all (virtual) orbitals of the previously gen- 
erated layers. These substitutions were augmented by a 
small subset of dominant single and double substitutions 
from core and valence shells, with the further restric- 
tion that at most one electron may be promoted from 
the core shells (i.e. in a double substitution at least one 
electron is promoted from a valence shell). All configu- 
rations from earlier steps were retained, with all previ- 
ously generated orbitals fixed, and all new orbitals made 
orthogonal to all others of the same symmetry. The ini- 
tial shapes of radial orbitals were obtained by means of 
a Thomas-Fermi potential, and then driven to conver- 
gence with the self-consistency threshold set to 10 -10 for 
spectroscopic orbitals and 10~ 8 for virtual orbitals, re- 
spectively. All radial orbitals were separately optimized 
for each atomic state. The Optimal Level form of the 
variational expression [32| was applied in all variational 
calculations. 



B. Configuration-interaction calculations 

In the third phase of the computations, the confi- 
guration-interaction calculations (i.e. without changing 
the radial shapes of the one-electron spin-orbitals) were 
performed, with multiconfiguration expansions tailored 
in such a way, as to capture the dominant electron corre- 
lation contributions to the expectation values. All single 
and double substitutions were allowed from several core 
shells and from both valence subshclls (i.e. 7s7p, or 7s6d, 
depending on the state) to all virtual shells, with the 
same restriction as above, i.e. that at most one electron 
may be promoted from core shells. The virtual set was 



reduced matrix element 
3 Pi VVW \\ 3 D2\ from 



TABLE I: The values of the 
of the Schiff operator 4\/3 7r 

Eq. (|10|l [a.u.]. Electron substitutions from different sets 
of spectroscopic orbitals: 6p[7s6d|7s7p] (second column), 
6sp[7s6d\7s7p] (third column), ... , 5spd6sp[7s6d\7s7p] (sixth 
column). Electron substitutions to different sets of virtual 
orbitals: (lv) one layer of virtual orbitals, (2v) two layers of 
virtual orbitals, ... , (5v) five layers of virtual orbitals. 'DF' 
= uncorrelated Dirac-Fock value. 



DF 


674 


virtual 




shells opened for substitutions 


set 


6p 


6sp 


5d6sp 


5pd6sp 


5spd6sp 


lv 


1931 


1892 


1862 


1749 


1647 


2v 


4597 


7544 


7553 


7598 


7754 


3v 


5141 


8073 


8112 


8125 


8327 


4v 


4760 


7936 


7947 


7889 


8117 


5v 


4642 


7834 


7818 


7743 


8217 



systematically increased from one to five layers. In a sim- 
ilar manner, several core subshells were systematically 
opened for electron substitutions — from the outermost 
6p up to the 5s5p5d6s6p subshells. 

The convergence of the calculations can be observed 
by monitoring the dependence of the matrix elements on 
the size of the virtual set as well as on the number of core 
subshells that are opened for electron substitutions. The 
effects of substitutions from AsApAdif and still deeper 
shells were estimated in our previous papers and turned 
out to be below 1 percent in the case of hyperfine struc- 
tures [42|, HH and a fraction of a percent in the case of 
transition rates [4J, W5| . 



IV. RESULTS 



Table[T] shows the values of the reduced matrix element 



3 Pi V 1 p (r) 



3 D, 



from 



of the Schiff operators 4V3 n 

Eq. (fT0)l . The expectation values were calculated for the 
multiconfiguration expansions obtained from all possible 
combinations of virtual sets, and from opening sequen- 
tially the core subshells, as described in section UlI Bl 

Table |TT] presents the values of the reduced matrix el- 
ement of the EDM operator V3 [ 3 Pi ||— er\\ 3 Z?2] from 
Eq. As in Table [H the EDM matrix elements were 
calculated for the multiconfiguration expansions result- 
ing from all possible combinations of virtual sets and from 
opening up the core subshclls, cf. section UlI Bl 

The data from Tables [T] and |TT] are collected also in 
Figure [H where they are presented as functions of the 
size NCF (see Eq. [1]) of the multiconfiguration expan- 
sion. Figure 2(a) shows the reduced matrix element of 
the Schiff operator, while Figure |2(b)| shows the value of 
reduced matrix element of EDM operator. The resultant 
EDM, induced by the nuclear Schiff moment in the 3 £>2 
state of the isotope Ra-223 is presented in Table HTT1 
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size of MCDHF space size of MCDHF space 

(a) Schiff (b) EDM 

FIG. 2: The values of the reduced matrix elements of the (a) Schiff and (b) EDM operators as functions of the size of the 
multiconfiguration expansion 



TABLE II: The values of the reduced matrix element of the 
EDM operator Vs [ 3 Pi ||-er|| 3 D 2 ] from Eq. © [a.u.]. The 
notations are the same as in Table I. 



DF 


5.849 


virtual 




shells opened for substitutions 


set 


6p 


6sp 


5d6sp 


5pd6sp 


5spd6sp 


lv 


4.818 


4.811 


4.810 


4.810 


4.810 


2v 


4.514 


4.467 


4.460 


4.460 


4.460 


3v 


4.270 


4.192 


4.157 


4.155 


4.155 


4v 


4.251 


4.168 


4.114 


4.109 


4.108 


5v 


4.240 


4.153 


4.083 


4.076 


4.075 



A. Electron correlation effects 

The calculations appear to be converged within 
the multiconfiguration approximation employed in the 
present paper. A comparison of the 'DF' value with 
the final result in Table IIIII shows, that correlation ef- 
fects are dominant in these calculations. An inspection 
of Tables Q] and HD reveals that the Schiff moment matrix 



TABLE III: The atomic EDM D A [W 9 S/I] of the \F,M F ) = 
1 3/2, 3/2} hyperfme state of the 3 I?2 level, induced by the 
nuclear Schiff moment in the isotope 2 HRa (/ = |) 
The notations are the same as in Table I. 



|) [a-u.]. 



DF 


0.0512 


virtual 




shells opened for substitutions 


set 


6p 


6sp 


5d6sp 


5pd6sp 


5spd6sp 


lv 


0.1812 


0.1181 


0.1163 


0.1092 


0.1028 


2v 


0.2694 


0.4378 


0.4374 


0.4399 


0.4489 


3v 


0.2850 


0.4393 


0.4378 


0.4383 


0.4492 


4v 


0.2627 


0.4293 


0.4244 


0.4208 


0.4329 


5v 


0.2555 


0.4224 


0.4144 


0.4097 


0.4346 



element is mainly responsible for the correlation correc- 
tion to the atomic EDM. Overall, the correlated Schiff 
value is more than an order of magnitude larger than 
the uncorrelated one. On the other hand, the correlated 
value of the EDM matrix element is about 30% smaller 
than the uncorrelated one. Together, they result in the 
atomic EDM, which is 8.5 times larger than the uncorre- 
lated value. There are several factors, which have to be 
taken into consideration in these calculations in order to 
capture the bulk of the electron correlation effects: 

1. Schiff interaction is localised inside nucleus 
(see Eq. [T4|) . therefore only the (one-electron) 



matrix elements 
■S1/2II V 1 p(r) \\p 3/2 



S1/2II V 1 p(r) \\pi/2 and 

contribute appreciably 

in Eq. (sec also" Table II in ref. [13). All other 
spin-orbitals may contribute through the indirect 
electron correlation effects (46j . 

2. Only the inner regions of the spin-orbitals con- 
tribute to the Schiff moment matrix elements. 
Therefore it is essential to accurately represent core 
polarisation in the atomic wave functions. 

3. For the electric dipole interaction, the outer parts 
of the spin-orbitals are important. Therefore the 
valence correlation effects have to be accounted 
for. The Babushkin and Coulomb gauge values dif- 
fer by 5 orders of magnitude in the calculation of 
the very weak 3 Pi — 3 I?2 transition [45[. This re- 
sults from the fact, that the transition energy is 
very small. Transition energy is difficult to repro- 
duce accurately in variational calculations, when 
the wave functions for the two levels are generated 
separately and when their one-electron orbitals are 
optimised independently. It is necessary to obtain a 
balanced description of both states |47fl , which can 
only be achieved if all two-body (double substitu- 
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tions) and possibly also three- and four-body (triple 
and quadruple substitutions) correlation effects are 
fully taken into account. In general, the two gauges 
exhibit different energy dependence: the Coulomb 
gauge is strongly energy dependent, therefore the 
Babushkin gauge value is usually adopted in un- 
saturated calculations, since it is less dependent 
on calculated transition energy value. The calcu- 
lations of energy level differences require well bal- 
anced orbital sets and typically demand highly ex- 
tensive multiconfiguration expansions. The results 
are usually in better agreement with experiment if 
a common set of orbitals is used for both states. 

4. On the other hand, the wave functions for the two 
levels have to be generated separately, if the ef- 
fects of non-orthogonality are to be taken into ac- 
count correctly Non-orthogonality between the 
spin-orbitals is essential in the evaluation of non- 
diagonal matrix elements (see Ref. [45| . 



B. Breit and QED corrections 

The calculations descibed in the previous sections as 
well as the data in Tables [IJ [Til and IIHI were obtained 
with only the Coulomb interaction included in the differ- 
ential equations which are iteratively solved during the 
self-consistent field process (SCF). The effects of mag- 
netic and retardation corrections can be evaluated by in- 
troducing the full Breit operator in the self-consistent 
field process. Such modified differential equations lead 
to change of the final wavefunction shape, which in turn 
modifies the matrix elements evaluated in the calculation 
of the Schiff moment. Quantum electrodynamics (QED) 
correction is dominated by two contributions. The self- 
energy (SE) part cannot be easily evaluated. Its con- 
tribution to the Csparity violation amplitude has been 
calculated in Ref. [48| and turned out to be of the order 
of -0.7 %. On the other hand, the vacuum polarization 
can be easily calculated by introducing the Uclhing po- 
tential into the SCF, as was done, for example, in calcu- 
lations of Li-like ions hyperfine matrix elements (49| or 
in Cs parity- violation amplitude 50]. In the latter case, 
it leads to a 0.4 % increase of the amplitude. If both 
are evaluated, the SE and VP contributions partly com- 
pensate. We have used the mdfgme code [13] in its 2008 
version, which now includes a Schiff moment option, to 
evaluate the effect of the Breit interaction and vacuum 
polarization on the Schiff moment. The results are pre- 
sented in Table QV] The inclusion of the Breit interaction 
in the SCF reduced the value of the Schiff moment by a 
factor of 1.7 %, while the vacuum polarization increased 
it by roughly 0.6 %. The total correction is a reduction 
of the order of 1.1 %. From Ref. [48[ one should expect a 
partial, mutual compensation of the VP and SE contri- 
butions. The total correction would then be dominated 
by the Breit contribution. 

The results presented in Table IIVI were obtained in 





Schiff 


Stark 


EDM 


correction 


DF 


674 


5.849 


0.0512 




Breit SCF 


664 


5.835 


0.0503 


-1.7 % 


Breit SCF + VP SCF 


668 


5.837 


0.0507 


-1.1% 



TABLE IV: Effect of the Breit and vacuum polarization on 
the reduced matrix element of the Schiff operator (Schiff), 
reduced matrix element of the EDM operator (Stark), and 
the atomic EDM {D A [10 9 S//]) of the \F,M F ) = |3/2, 3/2) 
hyperfine state of the 3 i?2 level, induced by the nuclear Schiff 
moment in the isotope 2 glRa (/ = §) [a.u.]. VP — vacuum 
polarization. DF — Coulomb interaction only. 

single configuration approximation, i.e., without electron 
correlation effects. Our experience indicates, that a fully 
correlated calculation, i.e., with the multiconfiguration 
expansion |T|) described in section IIII B[ and with the 
self-consistent Breit interaction would lead to a larger 
(absolute value of) total correction. However, the in- 
crease of the number of extra integrals involved (roughly 
by 2 orders of magnitude) and substantial convergence 
difficulties render such a calculation virtually impossible 
with today's computers. 

The final correction from Table [TVl (i.e., the reduction 
by a factor of 0.989 %) has been carried over to the final 
summary of our results, which is presented in Table fVl 



C. Accuracy 

Although we have full control of the core polarisation 
effects, the overall accuracy of these results depends pri- 
marily on the electron-correlation effects which were not 
included in the calculations, i.e. unrestricted double sub- 
stitutions and triple substitutions. These unrestricted 
substitutions had to be omitted due to software and 
hardware limitations. A similar approximation, based on 
single and restricted double substitutions, was employed 
also in our previous papers on radium [42l . l43l . 144 . [45| . 
The accuracy of the present calculations can be indi- 
rectly inferred from a comparison of the previously cal- 
culated transition rates and hyperfine constants with ex- 
periment. Reference [43[ showed hyperfine constants cal- 
culated in an approach similar to that employed in the 
present paper. As discussed there, the overall accuracy of 
the calculated magnetic dipole constants was 6 %, while 
for the electric quadrupole constants the estimated ac- 
curacy was 3 %. The interaction responsible for the hy- 
perfine shifts takes place in the vicinity of the nucleus. 
The bulk of the hyperfine integral/matrix elements come 
from the first oscillation of the one-electron wave func- 
tion [5l], H^j. Therefore, the calculated hyperfine con- 
stants depend heavily on the inner regions of the one- 
electron wave functions. The interaction of the electronic 
cloud with the nuclear Schiff moment also takes place in 
the vicinity of the nucleus. In fact, the bulk of the Schiff 
integral/matrix element comes from the nuclear skin (see 
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Eq. [T4|) . Therefore we might expect that the accuracy of 
the calculated Schiff moment matrix elements is compa- 
rable to the accuracy of the calculated hypcrfine con- 
stants. In conclusion, we might expect the accuracy of 
the calculated Schiff moment matrix elements of the or- 
der of 10 % or better. 

The electric-dipole moment matrix ele- 
ment (see Eq. I15[) has similar radial dependence as 
the electric-dipole transition elements, therefore their 
accuracy might also be expected to be of the same 
order. The calculated rate for the strong transition 
■"Pi — ^So showed excellent agreement between the val- 
ues calculated in Babushkin and Coulomb gauges 44j, 
which is an indication (but not a proof) of convergence. 
The accuracy of the calculated rate for the 3 Pi — ^So 
transition was also reasonably good [45[ . Although there 
was a 20 % difference between the values calculated in 
the two gauges, the Babushkin gauge value fell within 
the experimental limit [53]. The gauge difference was 
somewhat larger (35 %) for the calculated rate for the 
3 D2 — electric quadrupole transition [44| which is 
the only significant decay channel of the metastable 
state 3 D2- However, the calculated rates of other 
weak transitions were less accurate. In particular, the 
calculated rate for the 3 Pi — 3 D2 transition showed very 
large gauge dependence, which seems to indicate that 
the core-core correlation effects (which were omitted 
in present calculations) are important (see e.g. [HI). 
These correlation contributions arc currently beyond 
the capacity of the computer resources available to us. 
We were unable to include, or even to estimate the 
contributions of the omitted electron correlation effects, 
i.e., the unrestricted double substitutions and triple sub- 
stitutions, to the atomic EDM. Such calculation would 
require a multiprocessor cluster in excess of a hundred 
processors [55| . In order to estimate the uncertainty 
of the electric dipole moment calculation, we had to 
resort to an upper limit resulting from the hierarchy of 
electron-electron interactions. Table [II] demonstrates 
the saturation of the core-valence correlation correction, 
which is the dominant electron correlation effect beyond 
the Dirac-Fock approximation. Therefore we might 
expect that the entire (core-valence) electron correlation 
contribution would be a rather conservative estimate of 
the uncertainty of the electric dipole moment calcula- 
tions. In conclusion, we may assume the accuracy of the 
calculated electric dipole moment matrix elements of the 
order of 30 %. 

Another contribution that affects the overall accuracy 
of the results arises from the nuclear-density dependence 
of radial matrix elements in Eq. (jT4j) , i.e. from the distri- 
bution of the nuclear charge, as discussed in section Hi Bl 
However, the error of the calculated value of atomic EDM 
is dominated by the electric dipole moment calculation. 
Therefore we assumed 30 % accuracy of the final calcu- 
lated value of the Schiff moment enhancement factor, as 
shown in Tabic fVl 

A more stringent accuracy assessment for the present 



TABLE V: The atomic EDM induced by the nuclear Schiff 
moment in the 3 L>2 electronic state for three isotopes of ra- 
dium: 2 ! 3 8 Ra (I = \,F = §), 2 |lRa (/ = |,P = §), and 



2 g|Ra (/ = \-,F — |), respectively. Results are shown as 
functions of the size of the virtual orbital set, including elec- 
tron substitutions from the spectroscopic 5spd6sp[7s6d\7s7p] 
orbitals to: (lv) one layer of virtual orbitals, (2v) two lay- 
ers of virtual orbitals, ... , (5v) five layers of virtual orbitals. 
'DF' = uncorrelated Dirac-Fock value. Our values include 
the 'Breit SCF + VP SCF' correction from Table ITvTT he 
RHF+CI results in the last line are quoted from Ref. [TtJ. 



EDM [KfSyi] [a.u.; 



layer 


2 g|Ra 


g8 ±ia 


225 Rn 


DF 

lv 

2v 

3v 

4v 

5v 


0.0159 
0.0320 
0.1393 
0.1394 
0.1344 
0.1349 


0.0507 
0.1017 
0.4441 
0.4444 
0.4283 
0.4300 


0.0158 
0.0318 
0.1387 
0.1388 
0.1338 
0.1343 


final 


0.13(4) 


0.43(14) 


0.13(4) 


RHF+CI 


0.094 


0.30 


0.094 



calculations would require more experimental data to 
compare with. Hereby, we would like to encourage the 
experimenters to provide them. In particular, hypcrfine 
constants for other levels, as well as lifetime and tran- 
sition rate measurements for weak transitions would be 
very valuable. 

The results presented in the second and in the last col- 
umn of Table [V] illustrate the isotope-dependence of the 
atomic EDM. The two isotopes in question, 2 g|Ra and 
2 ||Ra, have the same nuclear spin (/ = 1/2), but they 
have slightly different nuclear shapes (see Eq. ([3]) and the 
subsequent discussion in section III Bp . The abovemen- 
tioned isotope-dependence, of the order of 0.4 %, arises 
primarily through the derivative dp/ dr in Eq. (|14[) . The 
isotope-dependence of the atomic wave functions (HJ has 
been neglected in the present calculations. 



V. CONCLUSIONS 

Radium is well suited for EDM experiments because 
of its large nuclear charge Z and the two low-lying levels 
3 Pi and 3 I?2 of opposite parity, which are separated in 
energy by only 5 cm -1 « 2 x 10~ 5 a.u. The 3 P/2 level 
is metastable with a very long lifetime of the order of 
4 seconds [45[ and therefore suitable for laser-ion traps. 
For radium, moreover, there are several isotopes with 
mass ranging from 209 up to 229, of which 2 g|Ra and 
2 ||Ra have nuclear spin I = 1/2. For EDM experiments 
the 7=1/2 isotopes are preferable since these isotopes 
cannot be disturbed by higher order electromagnetic mo- 
ments. 

The paper presents systematic computations of the 
static EDM of atomic radium in the 3 P2 level, induced 



10 



by the nuclear Schiff moment. Tabic [V] shows the atomic 
EDM induced by the nuclear Schiff moment in the 3 D 2 
electronic state for three isotopes of radium: 
^Ra (I = l,F= §,/* = 0.6133), 
2 iRa (I = |,-F = §, \i = 0.2705), and 
22 |Ra (I = |, F = = -0.7338). 
The wave functions for the two levels, 3 Pi and 3 D 2 , were 
generated separately, in order to correctly reproduce the 
effects of non-orthogonality between one-electron spinor- 
bitals. We demonstrate that core-valence electron cor- 
relation, which is the dominant electron correlation ef- 
fect beyond the Dirac-Fock approximation, contributes 
almost 90% of the total EDM value. Our final value is 
about 30% larger than the RHF+CI result from Ref. [l7t | . 
The difference can be atributed to different methods em- 
ployed to account for corc-valcncc electron correlation 
effects. Neither of these two calculations included the 
core-core correlation effects, which we believe to be the 
dominant source of uncertainty in our calculated value of 
the Schiff moment enhancement factor. 
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